#! /bin/bash
#Written by Graeme Eagles.
#Remove any existing files.
rm IPA.dat IPB.dat IPC.dat covariances.dat grnnam_covariances.dat grnnam_confint2d.dat grnnam_confint3d.dat aic.dat aai.dat LATESTIC LATESTAI

awk '{printf "%02d",$4-1}' C:/GMT/2019/GRNEUR3/inversion/data/iter.dat > LATESTIC
awk '{printf "%02d",$4-1}' C:/GMT/2019/EURNAM/inversion/data/iter.dat > LATESTAI

LATESTIC=`cat LATESTIC`
LATESTAI=`cat LATESTAI`

#assemble input files for covariance_add
awk '{if (NR>1) print $0}' C:/GMT/2019/GRNEUR3/inversion/data/poles_$LATESTIC.dat > aic.dat
awk '{if (NR>1) print $0}' C:/GMT/2019/EURNAM/inversion/data/poles_$LATESTAI.dat > aai.dat
paste aic.dat  C:/GMT/2019/GRNEUR3/inversion/data/covars_$LATESTIC.dat > ../grneur_input_for_covmatadd.dat
paste aai.dat  C:/GMT/2019/EURNAM/inversion/data/covars_$LATESTAI.dat > ../eurnam_input_for_covmatadd.dat
 
for CHRON in C13 C20 C21 C22 C24 C25

#echo "Combining matrices for the "$LATESTIC"th India-Capricorn and the "$LATESTAI"th Australia-India rotations."

do
echo "Working for chron "$CHRON "...."
#prepare file for rotation A
awk '{if ($4==chron) print $1,$2,$3,$5,$6,$7,$8,$9,$10,$11,$12,$13}' chron="$CHRON" ../grneur_input_for_covmatadd.dat > IPA.dat

#prepare file for rotation B
awk '{if ($4==chron) print $1,$2,$3,$5,$6,$7,$8,$9,$10,$11,$12,$13}' chron="$CHRON" ../eurnam_input_for_covmatadd.dat > IPB.dat

#prepare file for rotation C
awk '{if ($4==chron) print $2,$1,$3}' chron="$CHRON" ../grnnam_pol_circ.dat > IPC.dat

./covariance_add IPA.dat IPB.dat IPC.dat

awk '{print $0,chron}' chron="$CHRON"  covariances.dat >> grnnam_covariances.dat
awk '{print $0,chron}' chron="$CHRON"  confint.dat >> grnnam_confint2d.dat
awk '{print $0,chron}' chron="$CHRON"  confints.dat >> grnnam_confint3d.dat

done
